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A COMPARISON OF MEASURED AND CALCULATED UPWELLINC RADIANCE 
OVER WATER AS A FUNCTION OF SENSOR ALTITUDE 


Thom A. Coney, and Jack A. Salzman 

NASA- Lewis Research Center 
Cleveland, Ohio 


SUMMARY 

A comparison is made between remote sensing data 
measured over water at altitudes ranging from 30 m to 
15.2 km and data calculated for corresponding altitudes 
using surface measurements and an atmospheric radia- 
:ive transfer model. Data were acquired on June 22, 
1978 in Lake Erie, a cloudless, calm, near haze free 
day. Suspended solids and chlorophyll concentrations 
were 0.59i0.02 mg/1 and 2.42:0.03 pg/1 respectively 
throughout the duration of the experiment. Remote 
sensor data were acquired by two multispectral scan- 
ners each having 10 bands between 410 nm and 1040 nm. 
Calculated and measured nadir radiances for altitudes 
of 152 m and 12.5 km agree to within 167. and 147. re- 
spectively. The variation in measured radiance with 
look angle was poorly simulated by the model. It 
was concluded that an accurate assessment of the 
source of error will require the inclusion in the 
analvsis of the contributions made by the sea state 
and specular sky reflectance. 








INTRODUCTION 

Although various water quality parameters have been successfully measured 
using remote sensing techniques, there are a number of problems which must be 
solved before such techniques can be employed in operational monitoring. Theo- 
retical and laboratory studies, as well as field experiments, have demonstrated 
that concentrations of chlorophyll, suspended solids and other substances in the 
water can be interpreted from measurements of the intrinsic water color (refs. 1 
and 2) . While there is still work needed to quantify and expand these inter- 
pretative capabilities, of equal necessity is the development of techniques to 
accurately measure the water color from remote platforms. Aperture radiance 
values measured by remote sensors contain the necessary surface water color data 
but perturbations due to atmospheric interactions make their accurate extraction 
quite complex. Having the capability to correct for these atmospheric effects 
is essential before the unique perspective of satellite and aircraft systems and 
all thi-'lr ianerent advantages can be fully utilized. 

Several theoretical models and computer programs have been developed to per- 
form atmospheric corrections on remotely sensed data (refs. 3 and 4) but a major 
factor limiting their refinement and verification has been the lack of a suit- 
able data base with which to compare model calculations. During the summer of 
1978, the NASA Lewis Research Center organized a field experiment to acquire 
such data. These data have subsequently been used to test and refine an avail- 
able atmospheric radiative transfer model (ref. 5) in order to evaluate its pre- 
dictive capabilities. Other participants in this expei iment included the U.S. 


Environmental Protection Agency and Scripps Institute of Oceanography. Radiance 
data were acquired at altitudes ranging from 30 m (100 ft) to 15.2 km (50,000 ft) 
using two aircraft/multispectral scanner remote sensing systems. Optical data 
were also acquired both directly above the water's surface and at several sub- 
surface depths using a submersible spectroradiometer (owned and operated by 
Scripps). These data included the total downwelling irradiance and the upwell- 
ing radiance at three depths and the total downwelling irradiance in air. All 
measurements were made over a spectral range from 410 nm to 720 nm. 

The atmospheric radiative transfer model used calculates the total down- 
welling irradiance, the transmittance and the path radiance as a function of al- 
titude. The model includes the effects of multiple scattering by gases and par- 
ticulates in the atmosphere and of absorption by ozone in the near ultraviolet 
and visible part of the spectrum. Measured inputs to the model include date, 
time, barometric pressure, site latitude and longitude, aircraft altitude, and 
atmospheric optical thickness. Estimated inputs include the aerosol phase func- 
tion, ozone absorption, and standard aerosol optical depths above 5 km. 

Model predictions were compared with data acquired both at the surface and 
over the full range of altitudes. Downwelling irradiance was compared direct- 
ly with the surface measurements. The calculated upwelling radiances at the 
surface were compared with the aircraft scanner data by including model pre- 
dictions of atmospheric attenuation, path radiance and sky reflectance. 

THEORY 

In the formulation v'red here, the basic equation of remote sensing is: 

L(t.X,u 0 $ 0 ) - L 0 (t 0 .X.U 0 .« 0 ) T(t , A , u Q ) + L p (T.A,u 0 ,< 0 ) (1) 

uhere L(t , A , Uq , <»o) is the aperture radiance measured by the sensor, L 0 (x,A,ii 0 , 
$o) is the surface radiance, T(i,A,u 0 ) is the transmittance of the atmosphere 
and Lp(x ,A , m 0 ,^ 0 ) is the atmospheric path radiance. t,A,u 0 and $ 0 are the 
optical depth, wavelength, cosine of the solar zenith and solar azimuth angle 
respectively. T(t,A,p 0 ) and L p (t , A , p 0 , $ Q ) are determined by the atmospheric 
radiative transfer model. 

The radiative transfer model considered was developed by Dr. Robert Turner 
(ref. 5) under contract to Lewis Research Center. The model calculates the 
total downwelling irradiance, the transmittance and the path radiance for a 
plane-parallel, horizontally spatially homogenous, isotropic atmosphere as a 
function of altitude. The model includes the effects of multiple scattering by 
gases and particulates in the atmosphere and of absorption by ozone in the near 
ultraviolet and visible part of the spectrum. It is presently applicable to 
haze levels ranging from near fog conditions to the clearest possible atmosphere 
and the wavelengths between 400 nm and *A*2.2 m (with the exception of the strong 
absorption regions in the near IR) . 

The input variables can be separated into two categories, i.e., those that 
are estimated and those that are measured (see table I). Estimated values are 
used where a measurement would be difficult or where a variable is relatively 
constant. The model input format is such that a measured value can be used in 
place of an estimated value. 


EXPERIMENT 

Instrumentation 

The remote sensing data used in this experiment was obtained with two iaulti- 
spectral scanners, a ten band ocean color scanner (OCS) and an eleven band mod- 
ular multispectral scanner (M2S) . The ten bands on the OCS are between 425 and 
795 nm and have 20 nm bandwidths. Its full and instantaneous fields of view are 
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90° and 0.258° respectively. Its design altitude is 12.5 km. The OCS is mount- 
ed in a wing pod on the General Dynamics F- 106 shown in figure 1. 

The first eight bands of M^S are in the visible spectrum between A10 and 
720 nm and have bandwidths of about AO nm. Bands nine and ten are centered in 
the near infrared at 815 and 10A0 nm and have bandwidths of 90 nm. The eleventh 
band is in the thermal IR at 11 u. The full and instantaneous fields of view of 
the M*S are 90° and 0. 1A3° respectively. Its variable scan rate allow con- 
tiguous coverage at altitudes from 0.3 to 6.1 krr (1,000 to 20,000 ft). The M^S 
is mounted on the Convair C-131 shown in figure 2. 

The ship shown in figure 3 is the Environmental Protection Agency's R/V 
Crockett. On-board facilities include a wet laboratory and a chemistry labora- 
tory. This ship was used to obtain all surface truth data i.e., chlorophyll a, 
b, and c, phaeophytin, and suspended solids concentrations (ref. 6) as well as 
the optical properties of the water and the atmosphere. Using their submersible 
spectroradiometer (fig. A) Scripps Institute of Oceanography measured the total 
surface dowr.weiling irradiance, the subsurface total downwelling irradiance at 
three depths and the subsurface upwelling radiance at three depths; and used 
these data to calculate the surface upwelling radiance and the surface reflect- 
ance. The seven band solar radiometer shown in figure 5 was used to measure 
the optical thickness of the atmosphere. 

Experiment Procedure 

To evaluate the utility of an atmospheric scattering model as applied to re- 
motely sensed data it is most desirable to conduct the experiment under the 
least complicating environmental conditions; a cloud free atmosphere with a con- 
stant aerosol content, calm seas, and a large water area of low and unchanging 
chlorophyll and suspended solids concentrations. These conditions were most 
nearly attained on June 22, 1978. The scene shown in figure 6 was recorded at 
1A5A hours at the experiment site 1A km off Cleveland in Lake Erie. The sky 
was clear over the lake with only slight haze present and the seas were rela- 
tively calm. The suspended solids and chlorophyll concentrations were low, 
averaging 0.59±0.02 mg/1 and 2.A2±0.03 ug/1 throughout the day. 

Acquisition of Remote Sensing Data 

The direction and time of each aircraft flight line was determined by the 
orientation of the sun. Experience has shown that for solar zenith angles less 
than A0° the incidence of sun glitter from the wavelets increases significantly. 
To maximize available light however, data should not be obtained at zenith an- 
gles greater than 60°. These bounds established two windows on the 22nd of 
June (a morning window, 08A5 to 10A5 hours, and an afternoon window, 1625 to 
1825 hours) during which all remote sensing data were acquired. The C-131 ob- 
tained data at altitudes ranging from 30 m to 3.05 km (100 to 10,000 ft). A 
total of 21 ship overpasses were flown. The F-106 obtained data at altitudes 
ranging from 152 m to 15.2 km (500 to 50,000 ft) for a total of 10 overpasses. 

In all cases data were obtained with the aircraft flying in the solar plane and 
with the sun aft. 


Acquisition of Surface Truth Data 

The sequence of measurements made by Scripps began at 0915 and continued to 
1700 hours. The first measurement in the sequence was of the total downwelling 
irradiance made on deck with the submersible spectroradiometer (spectral range 
A10 to 710 nm) as shown in figure A. The spectroradiometer was then submerged 
under its support buoy (fig. 7) and the downwelling irradiance at three depths 
was measured. The spectroradiometer was then inverted and the upwelling ra- 
diance at three depths was measured. These subsurface measurements were then 
used to calculate surface upwelling radiance and surface reflectance (ref. 7). 
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The sequence was repeated such that a total of four measurements of the surface 
irradiance and eight measurements of the surface radiance and surface reflec- 
tance were made. Optical thickness measurements were made with two solar radio 
meters (fig. 5) every half hour from 0830 to 1800 hours. 

Water samples were obtained twice during each window and once before and 
once after solar noon (''-1330 hr) . Each time samples were obtained at three 
depths; 1 m below the surface, at the mid-Secchi depth and at the Secchi depth 

0 9 m) . 


RESULTS AND DISCUSSION 

A total of 31 overflights, 40 optical thickness measurements, 6 chlorophyll 
and suspended solids concentration measurements, and 26 irradiance and radiance 
measurements (using the Scripps submersible spectroradiometer) were made between 
0830 and 1800 hours. The purpose in acquiring this large amount of data was to 
estah’-sh the variability of each parameter throughout the day, to establish 
trends so that extrapolation in time was possible, and to acquire data at vary- 
ing solar zenith and azimuth angles. The following results are applicable to 
the general state of the water and air on June 22 and to thiee representative 
sensor altitudes i.e., 12.5 km (41,000 ft), the design altitude of OCS ; and 

3.05 km (10,000 ft) and 152 m (500 ft), the two altitudes coincident to both 
sensors . 

Using the Scripps measurements for L 0 in eq. 1, nadir aperture radiances 
L were calculated for comparison with the 30 and 152 m altitude remote sensor 
data. At these low altitudes T is nearly 1.0 and Lp is small relative to 
L 0 Therefore the effects of the atmosphere are considered minimal. Eleven com- 
parisons were made. In all cases the measured nadir radiances exceeded the 
calculated radiances (shaded area in fig. 8). The difference was attributed to 
specular reflectance of skylight at the air-water interface. The remote sensing 
equation was thus modified by the addition of a specular reflectance term L S T , 
i.e. , 

L - L 0 T + L p + L S T (2) 

where L s is here defined as p s E s (T,X,u 0 )/n. The specular reflectance p s is 
assumed to be independent of time and wavelength. E s (T,X,pp) is the diffuse 
sky irradiance and is calculated by the atmospheric scattering model. Figure 8 
identifies the three components making up the measured aperture radiance as a 
function of wavelength. These data were obtained at an altitude of 152 m at 
1736 hours. Note that at some wavelengths the specular radiance L s is com- 
parable in magnitude to the surface radiance L 0 . 

An estimate of the specular reflectance for June 22 was obtained by sub- 
stituting the radiances measured remotely from 30 and 152 m for L; the surface 
radiance measured by Scripps for L c ; and the calculated values for T, L p and 
E s into eq. 2 and solving for Pg. A value of 2.0±0.067. was obtained for the 
specular reflectance using a total of eleven low altitude data sets. This is 
consistent with accepted values for the specular reflectance. This value of the 
specular reflectance was used in all subsequent comparisons of measured and cal- 
culated radiance. An example of the calculated radiance for an altitude of 

12.5 km and 0915 hours and the relative magnitudes of the components of eq. 2 
are shown in fig. 9. Note the significant increase in the path radiance term. 

COMPARISON OF CALCULATED AND MEASURED RADIANCE 

The data from two representative altitudes and two solar zenith angles are 
presented in figures 10 and 11. The solid lines represent the measured nadir 
radiances and the symbols the calculated nadir radiances. The latter were 
based on the interpolated Scripps measurement of the upwelling radiance at the 
surface; on the transmittance, path radiance, and diffuse sky irradiance as de- 
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termined by the atmospheric scattering model; and on the estimated specular re- 
flectance. The error, averaged over all bands, was 14% and 167. at 12.5 km and 
152 m respectively. Since the error did not change significantly with an in- 
crease in altitude, it cannot necessarily be attributed solely to the model cal- 
culations. As presented in figures 10 and 11 there are four potential sources 
of error; the interpolation in time of the upwelling surface radiance as meas- 
ured by Scripps; the estimate of the specular reflectance; the calculation of 
the transmittance and path radiance by the atmospheric scattering model; and the 
accuracy of the radiance measurement at altitude. For this experiment the error 
associated with the last three is relatively fixed. However, greater confidence 
in the results can be realized if the comparison is made between surface re- 
flectances rather than between upwelling radiances. 

Because of the technique used by Scripps to measure the upwelling surface 
radiance, the measurement is most accurate when made near solar noon (ref. 7). 

All remote sensor data however were obtained at times greater than 2 1/2 hours 
from solar noon. To minimize the error due to interpolation of surface radiance 
it is possible to consider the surface reflectance rather than an absolute meas- 
ure of radiance. Surface reflectance is a function of water content but is 
nearly independent of time. Since the chlorophyll and suspended solids were con- 
stant throughout the day, the optimum single measure of surface reflect.. nee 
(also at solar noon) can be used independent of the time. (Note: the calcu- 

lated reflectance is strongly dependent on the model determination of the total 
downwelling irradiance. Figure 12 indicates the accuracy of this determination 
by comparing the calculated and measured irradiances for two times within the 
experiment windows.) 

Figure 13 shows the comparison between the surface reflectances calculated 
using the nadir value of the remote sensor data from three altitudes and that 
measured by Scripps at the water surface. The calculated surface reflectances 
are the averages from three overflights at each altitude. It is evident that 
the best agreement is realized for an altitude of 152 m where the effects of 
the atmosphere are small. And since the contribution from specular reflectance 
is significant at this altitude and less so at higher altitudes, the comparison 
indicates a need for some improvement in the correction for atmospheric scatter- 
ing. 


To determine the effect of the various input variables on the calculated 
surface reflectance, a parametric analysis was performed. Each variable was 
changed by some reasonable percentage from its best measured or best estimated 
value and used with the 12.5 km data. The resulting percent changes in the cal- 
culated surface reflectance are shown in table II. It should be noted that the 
magnitude of the parametric changes is very specific to the reference values 
about which the changes are made. A complete parametric analysis therefore 
should include additional reference conditions such as water with higher concen- 
trations of contaminants or an atmosphere having a greater level of haze. It is 
evident, however, that small changes in a combination of parameters could ac- 
count for the errors in the comparison shown in figure 13. 

Another dimension of the atmospheric scattering model that must be consider- 
ed, is its ability to correct for limb brightening, i.e., the increase in aper- 
ture radiance as j function of look angle. Figures 14 and 15 show the compari- 
son of the measured and calculated radiance as a function of look angle for an 
altitude of 152 m and 12.5 km. The calculated radiance was determined assuming 
the surface radiance and the specular reflectance to be independent of look an- 
gle. All of the calculated limb brightening is thus provided by the model. As 
can be seen the calculated limb brightening is much lower than measured. A 
parametric analysis indicated that no change in any one or a combination of in- 
put variables could account for this difference. To establish the cause of this 
difference will require a better understanding of the angular dependence of the 
upwelling surface radiance, the specular reflectance and/or the diffuse sky ra- 
diance . 
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SUMMARY OF RESULTS 


The purpose of this experiment was to assess the accuracy of atmospheric 
corrections to remotely sensed multispectral scanner data through use of the 
subject radiative transfer model. The conclusions resulting from this assess- 
ment are as follows: 

1. At an altitude of 152 m (500 ft) where the contribution from the atmo- 
sphere was not a factor, the surface measured radiance was consistently lower 
(about one-half) than the remotely measured radiance. The difference in the 
measured values was attributed to specular reflectance of skylight at the water 
surface . 

2. For the conditions reported here, the specular reflectance was determined 
to be 2% which is consistent with accepted values. Therefore, the calculated 
radiance at altitude was increased by an amount dependent on this percentage, 
the altitude and the calculated sky irradiance. The radiance thus calculated 
agreed to within 147. and 16% for sensor altitudes of 15.2 km (41,000 ft) and 

152 m (500 ft) respectively. 

3. The surface reflectance calculated using the remote sensor data measured 
at an altitude of 152 m compared well with that measured at the surface by 
Scripps . (Here surface reflectance is the ratio of the irradiance above the 
surface due to upwelling from below the water surface and the total downwelling 
irradiance above the surface.) However, the calculated reflectance increased 
as the sensor altitude increased. Small errors in a combination of atmospheric 
scattering model input variables could account for these increases. 

4. The remotely measured radiance compares poorly with the calculated ra- 
diance as a function of look angle. There is a strong indication that no change 
or combination of changes in the model input variable will account for this 
discrepancy . 

5. An accurate assessment of the corrections for atmospheric effects pro- 
vided by the radiative transfer model cannot be made without a better under- 
standing of the contributions made by the air-water interface, the sea state, 
and the sky radiance. Additionally the assessment should be made using data 
from both high and low concentrations of chlorophyll and/or suspended solids in 
the water. 
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Table I. Input Variables 


Estimated 

Extraterrestrial solar 
irradiance 

Barometric pressure 
at sensor altitude 

Background albedo 

Haze type and 

refractive index 

Ozone column 
density 

Ratio of aerosol scattering 
optical thickness to 
total aerosol optical 
thickness 

Optical depth above 
5 km 


Measured 

Date 

Time 

Experiment site 

latitude and longitude 

Sensor and site 
altitudes 

Surface barometric 
pressure 

Sensor center wavelengths 

Optical thickness as a 
function of wavelength 

Angle of the scan plane 
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Table II. Sample Parametric Analysis Using the 


Atmospheric Scattering Model 


Input 

variable 

Change in 
input variable, 
percent* 

Resulting change 
reflectance , 

425 nm 502 nm 

in calculated 
percent 

591 nm 672 nm 

Measured 
radiance 
at 12.5 km 

+10 

+81 

+48 

+21 

+65 

Skylight 

reflectance 

+50 

-12 

-12 

-6 

-35 

Time 

-6 min 

(from 0915 hr) 

+17 

+8 

+5 

+ 7 

Optical 

thickness 

+20 

-46 

-42 

-24 

-88 

Altitude 

-1 

+2 

+1 

+1 

0 

Background 

albedo 

-50 

+91 

+19 

+6 

+37 

Barometric 
pressure at 
altitude 

+20 

+13 

+4 

+1 

e 

Refractive 

index 

1.5 - 0.1 i 
(from Haze M 
1.5 - 0 i) 

+117 

+14 

+5 

+42 

Haze type 

Haze L 

(from Haze M) 

-4 

-4 

-2 

-9 

Difference in measured re- 
flectance for a sensor alti- 
tude of 12.5 km as compared 
with the reflectance measured 

+462 

+35 

+112 

+115 


at the surface, percent 
Unless otherwise noced. 
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Figure 1. F 106 high altitude aircraft - platform for Ocean Color 
Scanner. 
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Figure 2. C 131 medium altitude aircraft - platform for modular 
multispectral scanner. 




Figure 3. R/V Crockett - platform for surface truth 
Instrumentation. 


Figured Scripps Institute of Oceanography's submersible 
spectroradiometer. 
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Figure 5. Seven band solar radiometer. 



Figure 6. Typical sea and atmospheric state on June 22, 1978 
at experiment site. 
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30 meters from ship. (L^T). path radiance (Lp ). and surface radiance 

(LqT) for a sensor altitude of 152 m. 
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CALCULATED NADIR RADIANCE AT 12.5 km 
CALCULATED NADIR RADIANCE AT 152 m 
MEASURED RADIANCE AT NADIR 



WAVELENGTH, nm 

Figure 11. - Comparison of measured and calcu- 
lated nadir radiance ai 12.5 km and 152 m for a 
solar zenith angle at 
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Figure 15. - Comparison of measured and calculated 
radiance as a function of look angle for a sensor 
altitude of 12.5 km. a wavelength of 462 nm and a 
time of 0915 hr. 


